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Abstract. Surveys of weak gravitational lensing of distant galaxies will be one of the 
key cosmological probes in the future. We study the ability of such surveys to constrain 
neutrino masses and the equation of state parameter of the dark energy, focussing on 
how tomographic information can improve the sensitivity to these parameters. We also 
provide a detailed discussion of systematic effects pertinent to weak lensing surveys, 
and the possible degradation of sensitivity to cosmological parameters due to these 
effects. For future probes such as the Large Synoptic Survey Telescope survey, we 
find that, when combined with cosmic microwave background data from the Planck 
satellite, a sensitivity to neutrino masses of cr(^ mv) < 0-05 eV can be reached. This 
results is robust against variations in the running of the scalar spectral index, the 
time-dependence of dark energy equation of state, and/or the number of relativistic 
degrees of freedom. 
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1. Introduction 

With the advent of precision measurements of the cosmic microwave background (CMB), 
large scale structure (LSS) of galaxies, and distant type la supernovae (SNIa), a new 
paradigm of cosmology has been established. In this new standard model, the geometry 
is flat so that f2 to tai = 1, and the total energy density is made up of matter (Q m ~ 0.3) 
[comprising of baryons (Q b ~ 0.05) and cold dark matter (f2 CDM ~ 0.25)], and dark 
energy (f^e ~ 0.7, with equation of state w = Pj p ~ —1). With only a few free 
parameters this model provides an excellent fit to all current observations [1-4]. In 
turn, this allows for constraints on other, nonstandard cosmological parameters. 

Two very interesting possibilities discussed widely in the literature are: (i) a 
subdominant contribution to the total energy density in the form of neutrino hot dark 
matter (HDM) and hence limits on the neutrino mass m u , and (ii) an equation of state 
for the dark energy component that deviates from w — — 1 and/or is time-dependent. 

In the first case, data from atmospheric and solar neutrino oscillation experiments 
strongly suggest two important mass differences in the neutrino mass hierarchy, 
Am 2 tm ~ \Amj 3 \ ~ 2.4 x 10~ 3 eV 2 and Am 2 un ~ Am 2 12 ~ 7.9 x 10" 5 eV 2 (e.g., [5]). 
The simplest interpretation of this data set implies mi ~ 0, m 2 ~ \J Am^ and 

m3 ~ ^Am^m, i.e., the so-called normal hierarchy. Alternatively, since the sign of 
Am\ z cannot be fixed by atmospheric neutrino oscillations, an inverted hierarchy with 
m 3 ~ 0, m 2 ~ \J Am^tm and m x ~ \J Am 2 &tra is also allowed. Yet a third possibility is 

that the three mass eigenstates are nearly degenerate, with m x ~ m 2 ~ m 3 ^> ^Am^. 
If the last case prevails, then oscillation experiments can offer no useful information on 
the absolute neutrino mass scale. 

Cosmology provides an interesting measure of the neutrino mass through the 
neutrino's kinematic effects on large scale structure formation (e.g., [6]). Qualitatively, 
free streaming of a HDM particle causes essentially all of its own perturbations to be 
erased on scales below its free streaming length Af s . This erasure is simultaneously 
propagated to other species (eg., CDM) via the metric source term. The net effect 
on the present day large scale matter power spectrum P m (k) is a suppression of order 
AP m (k)/P m (k) ~ —8 Q u /Q m (assuming Q u <C Q m , and normalisation at k — > 0) at 
large wavenumbers k ^> k{ s ~ 27r/Af s , where Q u h 2 = yjm^/(92.5 eV). Indeed, within 
the ACDM framework, LSS observations to date have already been able to constrain 
the sum of the neutrino mass to J2 m u < 0.3 — > 1 eV (2a) [7-15], depending on the 
data set used. This limit should be compared with the current laboratory bound from 
tritium (3 decay experiments, m v < 2.2 eV (2cr) [16, 17], and the projected sensitivity of 
the upcoming KATRIN experiment, ~ 0.2 eV [18, 19]. 

In the second case, the simplest candidate for dark energy is a cosmological 
constant A, obeying w — — 1. A minimal extension of this is to assume that the 
dark energy equation of state is given by P = wp, where w is constant in time. Using 
this framework, various cosmological data sets have been analysed in many independent 
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studies (e.g., [20-22]). With the most recent data w is constrained to be extremely close 
to —1, the value for a pure cosmological constant (e.g., [23]). 

However, in many particle physics models of dark energy, the equation of state is 
expected to depend on the scale factor, and many authors have analysed the data in 
terms of dark energy models with varying equations of state (e.g., [24-26]). From an 
observational point of view, the simplest possibility is that w undergoes a transition 
between two asymptotic values. This is very often also the case in particle physics 
models of dark energy such as quintessence [27-29]. The same is also true in models 
where dark energy is a consequence of modified gravity [30,31]. Cosmological data to 
date show no evidence for time variations in the dark energy equation of state. However, 
for future high precision observations, it is natural to search for such time dependences. 
Later on when we analyse the ability of future cosmological probes to constrain neutrino 
masses and the dark energy equation of state, we will also include this possibility using 
a simple parameterisation. 

Recently, it was shown that presently available data admit an apparent degeneracy 
between m u and w [32]. Specifically, by extending w (taken to be constant in time) into 
the phantom regime (w < —1), the neutrino mass bound relaxes to J2 m u < 1-48 eV 
from a combined CMB+LSS+SNIa analysis, compared to J2 m u < 0.65 eV for the case 
of a fixed w — — 1 using the same data set. This m v ,w degeneracy can be explained as 
follows. The suppression of power at small scales due to a large m u can be accommodated 
by LSS data if Q m is simultaneously increased to "restore" the matter power spectrum to 
its canonical spectral shape. On the other hand, w and Qdc (and hence Q m = 1 — fide in a 
flat universe) are degenerate parameters in the Hubble expansion rate at redshifts probed 
by SNIa, such that an increased Vt m can be compensated for by a decreased w. These 
two degeneracies then combine to produce an apparent interplay between m u and w. 

It has been shown in [15] that the m u ,w degeneracy can be effectively broken by 
adding to the likelihood analysis new data from the Sloan Digital Sky Survey (SDSS) 
measurement of the baryon acoustic oscillation (BAO) peak [33], bringing the present 
bound on the sum of neutrino masses to the 0.3 — > 0.5 eV range. In this work, however, 
we consider the prospects for breaking this degeneracy with an alternative approach, 
namely, weak gravitational lensing. 

Weak gravitational lensing (or cosmic shear) of distant galaxies by the intervening 
large scale structure provides an unbiased way to map the matter distribution in the 
universe (e.g., [34-36]). Perturbations in the density field between the source and 
the observer induce distortions in the images of source galaxies. By measuring the 
angular correlation of these distortions, one can in principle reconstruct the underlying 
three-dimensional matter power spectrum, without making assumptions about the 
luminous versus dark matter bias (cf. conventional galaxy surveys). Since its first 
detection in 2000 [37-40], the sensitivity of weak lensing surveys has improved to 
the stage where they are able to provide constraints on the matter density Q m and 
the fluctuation amplitude cr 8 that are competitive with other cosmological probes 
(e.g., [41,42]). Several larger surveys are now in operation or in planning. These 
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include the Canada-France-Hawaii Telescope Legacy Survey (CFHTLS) [43], the Dark 
Energy Survey (DES) [44], Pan-STARRS [45], and VISTA [46], culminating in the 
most ambitious SuperNova/Acceleration Probe (SNAP) [47] and Large Synoptic Survey 
Telescope (LSST) [48] surveys. 

Furthermore, all future weak lensing surveys will provide photometric redshift 
information on the source galaxies. This additional information allows one to study 
tomographically the growth of the intervening large scale structure and the distance- 
redshift relation [49]. Indeed, weak lensing tomography is expected to reach a 
level of sensitivity comparable (and orthogonal) to SNIa for dark energy parameter 
determination (e.g., [50-56]). 

The structure of this paper is as follows. We outline the principle of weak lensing 
tomography and present the formalism for computing the convergence power spectrum in 
section 2. In section 3, we discuss some of the systematic effects relevant for tomographic 
studies. The Fisher matrix formalism and the projected errors for various cosmological 
parameters are presented in section 4. Section 5 contains our conclusions. 



2. Weak lensing tomography 

Distortions in the images of galaxies induced by density perturbations between the 
source and the observer can be described by the distortion tensor (e.g., [34,36]), 

dOi \ -72 -k + 7i / 

where 6 S and 6 denote the angular position of a light ray in the source and the 
image planes respectively, and the tensor ipij decomposes into a convergence (i.e., 
magnification, k < 1) and two shear (i.e., stretching, 71,72 <C 1) components. In 
the Born approximation, the convergence component is given by a weighted projection 
of the density fluctuations 5 = Sp/p along the line of sight, 

fXh 

k(0)= d X W( X )S(x,xO), (2.2) 
J 

where x — So dz/H(z) is the comoving radial distance, H(z) the Hubble parameter, 
and Xh denotes the distance to the horizon. The assumption of a flat universe is implicit. 
The lensing weighting function is defined as 

W( X ) = —a-\x)x r dx'nUx')^^ , (2-3) 

with Wo = (3/2) VL m Hl, and n ga i(x) = H(z) n ga \(z), where n ga \(z) is the number 
of source galaxies between redshift z and z + dz per steradian, and satisfies the 
normalisation condition / °° dz n g3 \(z) = n ga i. It is generally assumed that n ga i(z) takes 
the form 

' ' ' - ' ' (2.1) 







■(J exp 


[-©1 
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where the parameters a, (3 and z are chosen to fit the observed galaxy redshift 
distribution. Here, we fix a = 2 and (3 = 2. The choice of the parameter zq depends on 
the type of weak lensing survey: for a ground-based "wide" survey, z® — 1 corresponds 
to a median galaxy redshift z mcd ~ 1; for a space-based "deep" survey, the choice of 
zq = 1.5 gives z med ~ 1.5 (e.g., [57,58]). 

Where photometric redshift information on the source galaxies is available and the 
sources separable into tomography bins, equations (2.2) and (2.3) can be generalised to 
give the convergence of the ith bin via 

k(0) -»• «i(0), W(x) -> Wi(x), n g3l (z) -> m(z), n gal -> , (2.5) 

where n« = J °° cfe Tij(z) is the number density of source galaxies in the ith tomography 
bin, and J2i n i( z ) — n ga,i( z ) by construction. A simple model for rii(z) would be 
rii(z) = Q(z i+ i — z)Q(z — Zi)n ga \(z), where Zi and Zi+i denote the lower and upper 
limits of the ith tomography bin. 

2.1. Lensing observables 

Since lensing by large scale density perturbations leads to image distortions typically of 
order 1 % or less (much smaller than the intrinsic scatter in galaxy shapes), cosmological 
weak lensing must be measured statistically with a large ensemble of sources and lenses. 
By decomposing the convergence field into two-dimensional Fourier modes, 

«-(0) = ^^exp(i£-0), (2.6) 
i 

one can define the angular power spectrum of the convergence at multipole t between 
the ith and jth tomographic bin as 

<«j k?> = {^)H{t - £')Ci j . (2.7) 

In the Limber approximation (where only modes with longitudinal wavenumbers 
k 3 < x _1 < £/x contribute to the correlation), the convergence power spectrum is 
related to the three-dimensional matter power spectrum P m {k) via [59,60] 

rxh 

C 1 / = J o d X Wiix) Wjix) X' 2 Pm (k = £/x, z) , (2.8) 

assuming a flat universe. Thus weak lensing tomography with n t bins yields n t auto 
(i = j) and n t (n t — 1)/2 cross (i ^ j) correlation spectra. The case of n t — 1 corresponds 
to lensing without tomography, and we shall be using these two terms interchangeably 
throughout the text. 

Figure 1 shows an example of the convergence auto and cross power spectra 
expected for tomography with two bins. It is important to note that these power 
spectra are strongly correlated, and in their entirety do not provide as much information 
as a naive interpretation might suggest. The gain in information from tomographic 
binning can be quantified by the correlation coefficient between the power spectra of 
the individual tomography bins, R 1 ^ = C 1 / / \Jcf [49]. A strong correlation (i.e., 
~ 1) between two bins means that further subdivision is futile. 
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Figure 1. Convergence auto and cross power spectra for the standard ACDM 
model, assuming the galaxy distribution function (2.4) and two tomography bins at 
< z < 1.5 and 1.5 < z < 3. The la error bands are computed for a lensing survey 
with sky coverage / s k y = 0.01, ri ga j = 30 = 0.4, smeared over 

bands of width A£ = i/A. See section 4 for definitions. 



Note that while we have discussed only the convergence, existing and future weak 
lensing surveys will measure the source galaxies' ellipticities and hence the two shear 
components. However, by decomposing the shear signal into a gradient ("E'-mode") 
and a curl ("S-mode") component, it can be shown that the statistics of the .E-mode 
shear and of the convergence are identical (e.g., [61]). On the other hand, the S-mode is 
expected to yield a vanishing spectrum and therefore serves as a powerful test for residual 
systematics [62,63]. Henceforth, we shall consider the convergence power spectrum as 
representative of the lensing signatures. 



2.2. Dark energy and neutrinos 

The weak lensing signal (2.8) comprises of two distinct components: the lensing weights 
W(x) and the matter power spectrum P m (k). While the latter quantifies the universal 
matter distribution, the role of the former, in broad terms, is to project this three- 
dimensional distribution onto the two-dimensional sky for the observer. In general, for 
a single source at comoving distance Xo, the function W(x) peaks at roughly Xo/2. 
Thus, within the ACDM framework, a typical source distribution such as (2.4) peaking 
at z ~ 1 will be most efficiently lensed by structures at z ~ 0.5 (cf. Figure 2). 

The effects of dark energy are manifest in both the lensing weights and the matter 
power spectrum. In the case of W(x), the dark energy density f2 de and equation of state 
w(z) determine the Hubble expansion rate at low redshifts, 

H 2 (z) = Hi [(1 - fi de ) (1 + zf + n de e 3xJ>'[i +w( *')]l ? (2.9) 
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and hence the distance-redshift relation, x — lo dz/H(z), responsible for the projection 
(cf. Figure 2). In the above equation we have assumed spatial flatness. This assumption 
will be kept throughout the analysis. 

On the other hand, dark energy suppresses the growth of structure relative to the 
case of a flat, matter-dominated universe. In the absence of massive neutrinos, the linear 
matter power spectrum evolves as P m (k,a) = D 2 (a)P m (k,a = 1), independently of k. 
Here, D(a) is the growth function normalised to D(a — 1) = 1, and can be obtained by 
solving the differential equation (e.g., [64]) 

2 7JT-^ + t 5 " 3w(a)fi de (a)]-^- + 3[1 - w(a)]Q de (a)g(a) = 0, (2.10) 
a in cr a ma 

where g(a) = D(a)/a, and fi de (a) is the dark energy density at epoch a. The effects 

of fide and w on the growth function are generally mild compared to those on W(x)- 

For the case of w — — 1.2, for example, the lensing weight W(x) increases by a uniform 

3^4% relative to the w — — 1 case across < z < 2, while the change in the growth 

function D(a) is a minute —0.05 % at z — 0.5 at the peak of the lensing efficiency and 

only reaches —2 % at z = 2 at the tail of W(x)- Thus the dominant constraints on the 

dark energy parameters will come from projection. In addition, fi de (= 1 — ^m) impacts 

immediately on the epoch of matter-radiation equality, which in turn determines the 

location of the turning point in the matter power spectrum. 

Neutrino masses in the sub-eV range do not play a role in the universal expansion 
rate at low redshifts. Their early free-streaming effect on the shape of the matter 
power spectrum is however well known. As mentioned earlier, the loss of power at 
large wavenumbers due to a large m v can be offset by simultaneously translating the 
matter power spectrum towards large values of k. This translation can be accomplished 
by increasing (decreasing) Q m (fide) and hence moving the matter-radiation equality 
epoch forward to an earlier time when the causal horizon is smaller. In this way, m v 
and fide form a degenerate pair. 

Lensing without tomography, like angular galaxy surveys, makes use of a projected 
signal accumulated over a large distance (~ 4000 h^Mpc). The information content 
of this signal, especially in the radial x direction, is therefore limited. Tomography 
binning allows one to recover some radial information through a set of lensing weights 
Wi(x) that peak at different values of x- This is particularly useful for establishing 
the distance-redshift relation and hence breaking the fi de ,w degeneracy. As shown in 
Figure 2, two sets of {fidc, - ^} with lensing weights that are indistinguishable to the 
naked eye in the case of no tomography begin to show perceivable signs of differences 
when three tomography bins are used. 

Tomography will not improve directly the neutrino mass bound. Even though 
massive neutrinos continue to affect the growth of small scale density fluctuations at 
low redshifts, the shape of the matter power spectrum evolves by a negligible amount 
over the z range probed by weak lensing. To illustrate this point, we note that for 
"Y^m v = 0.3 eV, the normalised power spectrum P m {k)/P m (k — > 0) changes by less than 
2 % at large values of k between z = and z = 1.5. This is smaller than the overall 
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X [h _1 Mpc] 

Figure 2. Top: Distance redshift relation for four sets of dark energy parameters 
{fide, w} listed in the Figure, where w is taken to be constant in time. Middle: Lensing 
weights W(x)/£l m , without tomography, for the same four sets of {fide, w}, assuming 
the galaxy distribution function (2.4) with a = 2, /? = 2, and zq = 1. Bottom: Lensing 
weights for three tomographic divisions at Q < z < 1, 1 < z < 2, and 2 < z < 3. 

power suppression expected even for a 0.07 eV neutrino today (~ 8 Q u /Q m ~ 5 %). 
Nonetheless, once the fl^e, w degeneracy is lifted by projection effects, the f2de, n^u shape 
degeneracy will be correspondingly reduced and the bound on m v tightened. 

3. Systematic effects 

A number of systematic effects have been identified that could severely limit the 
potential of weak lensing surveys for cosmological parameter extraction. We discuss 
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some of these below. For other systematics not covered here, we refer the reader to [36]. 



3.1. Photometric redshift uncertainties 

Future surveys will observe a projected billion galaxies, a sample too large for individual 
spectroscopic redshift measurements. By observing the galaxy broad band features (i.e., 
broad band photometry) and comparing them to predictions from galaxy spectral energy 
distributions, redshift information can be obtained much more quickly at the cost of 
lower accuracy. 

The uncertainties associated with the photometrically measured redshift are 
typically ±0.034 -> 0.1 (Hubble Deep Field [65], CFHTLS [42], SNAP with IR- 
filters [57], Combo-17 survey [66] with multi-colour filters etc.), compared with ±0.001 
using spectroscopic methods (DEEP2 [67], VVDS [68], etc.). Assigning redshifts to 
galaxies using photometric methods should not contribute significantly to the error 
budget in weak lensing (without tomography) measurements. For tomographic studies, 
however, binning with photometric redshifts z p ^ that suffer from uncontrolled scatter 
and bias relative to the true (spectroscopic) redshifts z can lead to considerable error 
degradation in the dark energy parameter extraction. The authors of [69, 70] found that 
the scatter and the bias per redshift interval 5z = 0.1 need to be known to an accuracy 
of 0.003 — > 0.01, in order not to degrade the errors in f2 dc and w by more than a factor 
of 1.5. 

We account for uncertainties in the photometric redshift following the approach of 
[69]. The true number of galaxy per steradian at redshift [z, z+dz] in the ith tomography 
bin is defined as 

m(z)= l dz ph n g3l (z) p(z ph \z) , (3.1) 
Jzt P h 

where the function p(z p i l \z) is the probability that a given z will be measured as 
z ph . This formulation assumes that n ga \{z) will be accurately determined by an 
independent, dedicated spectroscopic survey (an alternative formulation for the case 
of a photometrically derived n ga \(z) can be found in [70]). For simplicity, we choose 
p(zph\z) to be a Gaussian, 



p{z ph \z) = J_ exp 
V2%cr z 



(%>h ± ^bias Z) 

2^ 



(3.2) 



where cr z (z) and z hias (z) are the scatter and the bias parameters respectively. Following 
[69], we discretise o z (z) and z hias (z) into N pz = 31 sets of parameters {a*, zl ias }, where 
o l z = cr z (zi) and ^ ias = z hias ( z i), on an equally spaced grid from z = to z = 3. Later 
on, each of o~\ and z l Uas will be treated as a free parameter in our error forecast analysis. 

The exact forms of cr z (z) and z hias (z) depend on how the photometric redshift is 
measured by the lensing survey. For example, o z can be reduced by going to space (to 
eliminate the problem of absorption in the Earth's atmosphere), or by implementing 
IR-filters (sensitive to galaxies at 1 < z < 2). Another possibility is to calibrate the 
survey internally by doing spectroscopy on a fraction of the source galaxies [71]. 
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In this study, we choose as our fiducial model o~\ = 0.05 and ^ ias = 0, and for these 
values, we have verified that the accuracy requirements for the scatter and the bias are 
roughly consistent with those found in [69]. 

3.2. Shear measurement errors 

Technically, shear measurement is performed not on the actual galaxy images, but 
on images convolved with a point spread function (PSF) which describes the various 
image distortions originating from the telescope, the camera, and, in the case of 
ground-based surveys, the atmosphere. Errors in the lensing signal can therefore 
arise from uncorrected spatial and temporal variations in the PSF, as well as from 
the deconvolution and shear estimation method itself. Shear measurement errors can 
generally be discussed under the umbrellas of multiplicative and additive errors. 

3.2.1. Multiplicative errors A number of algorithms have been developed so far to 
estimate the shear in an observed image, all of which must be calibrated against 
simulations (see, e.g., [72,73] for a summary). Uncertainties in the calibration can 
lead to the measured shear being over- or under-estimated by a constant multiplicative 
factor, and hence impact directly on the Q m and erg parameter determination. Shear 
calibration errors are estimated to be 2% currently, and are expected to be controlled 
to 1% by the time SNAP and/or LSST start [74]. 

To incorporate this effect into our analysis, we follow the approach of [70], and 
define the theoretical convergence power spectrum as 



where is the irreducible part of the shear calibration error in the ith tomography bin, 
after averaging over all directions and redshifts within the bin. Thus, for tomography 
with n t bins, there are n t shear bias parameters. 

3.2.2. Additive errors Additive errors arise when artificial shear signals are induced by, 
e.g., anisotropies in the PSF. In general, the PSF varies both spatially and temporally, 
and the anisotropies incurred by these variations must be calibrated, usually using stars 
in the observed fields, for each individual image. Current simulations suggest that 
additive PSF errors can be controlled to under 10 -3 [74]. 

The effects of additive errors are more difficult to parameterise, since the errors 
will generally depend on I. A first attempt in [70] defines the total convergence power 
spectrum as C\ 3 = C\ 3 + C| dd , where the additive error spectrum C^ d is assumed to 
be a power law in t. Within this framework, the authors of [70] conclude that additive 
errors must be controlled to ~ 10~ 5 . We shall not be considering additive errors in 
this work, but simply note here that the size of current simulations are too small to 
determine if the 10~ 5 accuracy requirement is realistic [74]. 




(3.3) 
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Figure 3. Nonlinear matter power spectra constructed from the halo model for various 
neutrino masses (solid lines). Also shown are the corresponding power spectra from 
the linear theory (dotted lines). 

3. 3. Matter power spectrum 

The sensitivity of weak lensing surveys peaks at scales of several arcminutes, i.e., 
t ~ 1000 [75] (cf. Figure 1). At these scales the shear statistics is dominated by density 
perturbations of comoving wavenumbers k ~ 1 — ► 10 /iMpc" 1 at z ~ 0.5. Thus, weak 
lensing probes a convolution of linear and nonlinear scales at present and earlier times, 
and the linear power spectra P L (k) computed by, e.g., CMBFAST [76] alone will not 
suffice to determine total matter power spectrum P m (k). 

Assuming that baryons trace CDM exactly, accurate predictions of P m {k) can be 
obtained by means of collisionless iV-body simulations. At present various cosmological 
simulation codes agree at the 5 —>■ 10 % level [77] for ACDM cosmologies. In the 
future, however, P m (k) will have to be calibrated to about 1 — > 2 % accuracy in order 
not to compromise the cosmological parameter extraction from ambitious surveys such 
as SNAP and LSST [78]. Inclusion of baryon physics introduces an extra layer of 
complexity; dissipative processes such as radiative cooling and star formation enters at 
k > 1 /iMpc, and require in addition the use of hydrodynamical simulations to accurately 
model the behaviour of the baryon gas. We note in passing that large volume simulations 
of dark energy cosmologies to date have focussed mainly on ACDM, although there have 
been some efforts in the direction of to ^ -1 and/or a time- varying equation of state 
(e.g., [79-82]). To our knowledge, simulations of dark energy cosmologies with mixed 
dark matter have not yet been reported in the literature. 

A number of semi-analytical methods and fitting formulae have been devised [83- 
87], some of which tested to roughly 10 % accuracy against A^-body simulations [53, 88], 
to map between the linear and nonlinear power spectra in CDM cosmologies. Virtually 
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all weak lensing parameter constraint analyses to date employ one of these mapping 
methods. In this paper, we correct for nonlinearity using the halo model [84-86], 
assuming that it can be generalised to general dark energy and mixed dark matter 
cosmologies. Here, the total matter power spectrum is taken to be a weighted average 
of the neutrino (is), and the combined baryon (6) and CDM (c) power spectra, 



with weights fi = Qi/Q m , i = v,b,c, and Vt rn = Q u + Qf, + fi c - The nonlinear baryon 
and CDM power spectrum P^ c (k) is mapped from its linear counterpart, 



following the halo model prescription. See Appendix A for details. 

Figure 3 shows the nonlinear matter power spectra constructed in this fashion 
for various neutrino masses. Clearly, the suppression of power due to neutrino free- 
streaming extends well into the nonlinear regime (k > 0.2 /iMpc). Indeed, the halo 
model predicts an even stronger suppression on nonlinear scales than in the linear 
region (k < 0.2 hMpc). However, compared with their corresponding uncorrected power 
spectra (i.e., power spectra calculated from the linear theory, extended into the nonlinear 
regime), this relative suppression, i.e., the ratio AP m (k) / P m (k) , in the corrected spectra 
seems to be somewhat weaker. (Observe in Figure 3 that the gaps between the solid 
lines seem to be narrower than those between the dotted lines at a given k.) This result 
remains to be verified by simulations. 

We do not consider nonlinearity in the neutrino power spectrum P„{k) incurred 
from late-time infall into existing halos [89], since it affects the total matter power 
spectrum and hence the weak lensing observables only at the ~ 0.1% level [90,91]. 
We ignore also corrections due to cooling baryons and the intra-cluster gas [92,93], 
although results from recent iV-body/hydrodynamical simulations suggest that baryon 
physics affects the weak lensing signal at the < 0.5 % and the 1 — > 10 % levels at 
£ < 1000 and 1000 <£< 10000 respectively [94]. 

4. Error forecast 

We use the Fisher matrix formalism to estimate the constraints on various cosmological 
parameters from weak lensing tomography. 

4-1. Fisher matrix 

The Fisher matrix is defined as 



where L(x,p) is the likelihood function of the data set x given the parameters 
P — {Pi, ■ ■ ■ iPa) °f the theoretical model. When evaluated in the vicinity of the 




(3.4) 



Pt +C (k) = (fb + fc)- 1 hy/lfm + fc^PRk) 



(3.5) 




(4.1) 
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Table 1. Generic weak lensing survey parameters. 





/sky 


zo 


Trms 


rigai (arcmin 2 ) 


Wide, LSST-like 


0.7 


1.0 


0.4 


30 


Deep, SNAP-likc 


0.01 


1.5 


0.25 


100 



fiducial model, the Fisher matrix quantifies the best statistical error on the parameter 
p a achievable with a given data set via 

o-( Pa ) > v/C^U, ( 4 -2) 
where c(p a ) denotes the la error on p a after marginalisation over all other parameters. 

Neglecting non-Gaussian corrections, the Fisher matrix for weak lensing tomo- 
graphy has the form 

F«„ = pi + l/2)/, ky At T r (c^CT'g) • (4.3) 

*-min ^ ' 

Here, we take the lensing observable Ce to be the number weighted convergence power 
spectrum Ce = HifijC 1 / , since it contains more information than does the binned 
spectrum C 1 / to reflect changes in the photometric redshift parameters [69]. The 
variance Ce is correspondingly defined as a number weighted sum of the convergence 
and the noise power spectra, Ce = fiifij{C % ^ + hjlrms/^-i)- Note that in the case of fixed 
photometric redshift parameters, these definitions lead to the same F af3 as their more 
common representations, Ce = Cf and Ce = C\ 3 + o~ij7r ms /fii- We consider a variety of 
tomography bin numbers (n t = 1, 3, 5, 8), each with equal redshift binning Az between 
z = and z = 3 (e.g., Az = 0.6 for n t — 5). 

The minimum and maximum multipoles £ m i n and £ max are taken to be 50 and 3000 
respectively, the latter so as to avoid the highly non-Gaussian signals at the very small 
scales. Between £ min and £ max we use 100 logarithmic bins of width A£; variations in the 
£ binning have little effect on our results, since the convergence power spectrum is fairly 
featureless. The other parameters — the sky coverage / S k y , the surface density of galaxies 
^gai = Hi^ii and their rms intrinsic ellipticity 7 rms — are generally survey-dependent. In 
our analysis, we consider two generic weak lensing surveys: a ground-based LSST-like 
"wide" survey, and a space-based SNAP-like "deep" survey, with survey parameters 
displayed in Table I. 

We are primarily interested in how parameter constraints derived from CMB data 
can be improved with the inclusion of weak lensing. Since the primary CMB anisotropies 
and the galaxy shears probed by weak lensing surveys are generated at two vastly 
different epochs, we may take them to be uncorrelated. Thus the total Fisher matrix is 
simply the sum of two independent matrices, 

which is a good approximation provided we confine the analysis to the unlensed CMB 
anisotropies. (Weak lensing of the CMB and of galaxies are necessarily strongly 
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Table 2. Experimental specifications for Planck, assuming one year of observation. 



/sky 


v (GHz) 


9b (arcmin) 


A T (10" 6 T) 


A E (10- 6 T) 


0.65 


100 


9.5 


2.5 


4.0 




143 


7.1 


2.2 


4.2 




217 


5.0 


4.8 


9.8 



correlated since the lensing signals are generated by much of the same intervening large 
scale structure.) 

The Fisher matrix for the unlensed CMB anisotropics F™ B is given likewise by 
equation (4.3), but with Ce = Ce + N e , and 

_ / CJ T CJ E \ Ar / NJ T \ 

C ' = [<V* ch)> N <={ o aH' (45) 

where Cj T , Cf E , and Cj E are the auto and cross power spectra of the CMB temperature 
and E-type polarisation. These are computed using CMBFAST version 4.5.1 [76]. We 
do not consider 5-type polarisation, although it will most likely be measured by the 
next generation of CMB probes. The noise power spectra Np vary with the detector 
performance according to 



N% = (6 b A a ) 2 exp[£(£ + 1)02/8 In 2], (4.6) 



where v denotes the detector frequency channel, 0& is the beam width, and A a is the 
temperature/polarisation sensitivity per pixel. Table 2 shows the expected performance 
of the Planck satellite [95] (assuming one year of observation), and we take £ max = 1500 
in our analysis. 

Prior information on the photometric redshift uncertainties (i.e., the parameters 
z^ias an d cr*) and the shear calibration bias (i.e., /j) can be included in the analysis by 
adding the diagonal matrix F^ or = 5 a /3 <r(p a ) _2 to the total Fisher matrix (4.4). 

4-2. The fiducial model 

We consider a spatially flat cosmological model with eleven parameters, the fiducial 
values of which are as follows: J2 m u = 0.07 eV, Q dc = 0.7, Vt b h 2 = 0.0245, 
VL c h 2 = 0.1225, w = -1, w a = 0, r = 0.05, n s = 1, a s = 0, a 8 = 0.9, and N cS = 3. 
In this model, the time-dependent dark energy equation of state is parameterised as 
[96, 97] 

w(z) = w + w a [l -a(z)\. (4.7) 

See, e.g., [24-26] for alternative parameterisations. We also allow for nonstandard 
relativistic degrees of freedom N cS ^ 3, and a running scalar spectral index a s evaluated 
at ko = 0.05 WVlpc -1 . These extra parameters, {N e s, a s ,wo,w a }, are known to be 
degenerate with the neutrino mass m u . For simplicity, we take the number of massive 
neutrino species to be one. Spatial flatness implies Qdc + Q u + fib + Q c = 1- 
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Figure 4. Projected la contours in the dark energy and neutrino parameter space 
from an LSST-like ground-based wide lensing survey (with and without tomography) 
and Planck for an eleven parameter cosmological model. The red (dot-dash), blue 
(solid), green (dashed), and gray (dotted) lines correspond to one, three, five, and 
eight tomography bins respectively. Each contour comes from marginalising over the 
other nine cosmological parameters and all 2 x N pz + n t systematic parameters, with 
Gaussian priors imposed on the latter. 



In addition to the eleven cosmological parameters, we marginalise over 2 x N pz + n t 
nuisance parameters, cr*, zl ias , and /j. The fiducial values of these parameters are 
taken to be a\ = 0.05, z l Uas = 0, and fa = 0, and the Gaussian priors c{a l z ) = 0.01, 
cr(^ias) = 0.01, and cr(fa) = 0.02 are used in the marginalisation. We do not employ 
priors on the cosmological parameters. 
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Figure 5. Same as Figure 4, but from a SNAP-like space-based deep survey and 
Planck. 



4-3. Results 

Figures 4 and 5 show the projected la contours in the dark energy and neutrino 
parameter space expected from two generic weak lensing surveys, with and without 
tomography, in combination with prospective CMB data from the Planck mission. (See 
Table 3 for all projected errors, and Table 4 for the parameter error correlation matrix 
F a p/ \J F aa Fpi3.) Clearly, the addition of tomographic binning significantly reduces the 
errors for all four parameters, {J2 m u, w , w a , fide}- Focussing on the LSST-like wide 
survey (Figure 4), we see that the improvement on the w a error estimate is roughly a 
factor of five. For the neutrino mass, the gain is slightly more moderate: roughly a factor 
of three difference between the cases of no tomography and tomography with five bins. 
This follows from the fact that tomography does not directly add extra information on 
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Table 3. Projected la errors from various combinations of Planck and weak lensing. 
The labels "Wide" and "Deep" denote the two generic lensing surveys considered in 
this work, and the suffix number indicates the number of tomography bins used in the 
analysis. 





Planck 


+Wide-1 


+Wide-5 


+Dccp-1 


+Deep-5 


aQ2m v ) (cV) 


0.48 


0.15 


0.043 


0.39 


0.047 


a(£l dc ) 


0.08 


0.020 


0.0068 


0.036 


0.0099 


a(n b h 2 ) 


0.00028 


0.00016 


0.00013 


0.00024 


0.00014 


a(n c h 2 ) 


0.0026 


0.0017 


0.0015 


0.0019 


0.0015 


a(w ) 


0.83 


0.093 


0.034 


0.35 


0.045 




4.0 


0.39 


0.081 


1.7 


0.063 


a(r) 


0.0046 


0.0043 


0.0042 


0.0045 


0.0043 


a{n s ) 


0.0089 


0.0056 


0.0028 


0.0074 


0.0047 


a(a s ) 


0.024 


0.013 


0.0061 


0.020 


0.012 




0.084 


0.019 


0.0076 


0.030 


0.0092 




0.19 


0.11 


0.067 


0.14 


0.093 



the shape of the matter power spectrum (see discussion in section 2.2). 

The SNAP-like deep survey (Figure 5) offers even more significant relative 
improvement, since it probes further back in redshift where the effect of dark energy on 
the distance-redshift relation is more prominent (cf. Figure 2). However, owing to its 
limited sky coverage, the absolute errors (with tomography) are comparable to those 
from the wide survey. Should the value of / s k y be increased, one can scale the errors 
presented here roughly via 

0~(Pa, /sky) ~ (/sky.fid/ fskyY^^iPa, /sky.fid), (4.8) 

where / s k y ,fid denotes the sky coverage of the fiducial survey considered in this work. This 
expression is applicable for {J2 m„, Wo, w a , Qdc} with tomography, since the constraints 
on these parameters come predominantly from lensing, not from CMB. 

It is interesting to compare these constraints with those derived from CMB 
alone. Table 3 shows the projected errors on various cosmological parameters from 
Planck alone, and in combination with a selection of lensing surveys. CMB alone has 
very little sensitivity to the dark energy equation of state; lensing is crucial for the 
determination of w and w a . On the other hand, although the projected error on 
the neutrino mass from Planck alone is comparable to what can be achieved today 
with CMB+LSS+SNIa+BAO, weak lensing (without tomography) brings it down to a 
level comparable to the projected sensitivities of future galaxy surveys (e.g., [98]). The 
addition of tomography binning pushes the error down even further. 

Of course, these impressive sensitivities come at a price: appropriate priors have 
been imposed on the systematics parameters in order not to degrade the sensitivity 
to any one (cosmological) parameter by more than ~ 50 % relative to the best 
achievable error by a systematics-free lensing survey. Figure 6 shows the degradation 
in the projected sensitivities to {J] m„, w , w a , fide} as a function of priors on the shear 
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Figure 6. Degradation in the projected sensitivities, relative to the best achievable 
errors by a systematics-free lensing survey, as a function of priors on the shear 
calibration bias and photometric redshift errors. Contours of constant degradation 
are shown in intervals of 20 %, in ascending order from left/bottom to right/top. The 
thick black line in each plot indicates a sensitivity degradation of 30 %. Left. Contours 
for Planck+Wide-5. Right: Contours for Planck+Deep-5. 



calibration bias and photometric redshift errors for two combinations of Planck and 
lensing tomography with five bins. Clearly, the largely vertical contours suggest that 
the accuracy to which the photometric redshift uncertainties can be controlled is a 
crucial factor in deciding a survey's parameter sensitivities. For the neutrino mass 
determination, for example, the sensitivity degrades rapidly as the accuracy requirement 
for {cr* , ^biasl * s relaxed. In contrast, the degradation due to shear calibration errors 
saturates at ~ 70 % for Q^c (Planck+Wide-5), and much less for other parameters. 
This is qualitatively consistent with the findings of [70], where the shear calibration 
errors are said to enter a "self-calibration" regime. 

Changing £ max to a lower value degrades further the parameter sensitivities. For the 
neutrino mass extraction, the degradation is ~ 30 % from £ max = 3000 to £ max = 1000. 
For the dark energy parameters, on the other hand, the sensitivities degrade generally 
by less than 20 %. 

Focussing now on the m u ,w degeneracy, Table 5 compares various projected 
constraints on the neutrino mass for the eleven parameter model (cf. section 4.2) and 
a seven parameter model which assumes {N e s,a s ,wo,w a } to be fixed at {3,0,-1,0}. 
From Planck and weak lensing without tomography, the expected la error for the two 
models, 0.15 eV and 0.082 eV, differ by almost a factor of two. With the inclusion 
of five tomography bins, the errors reduce to 0.043 eV and 0.037 eV respectively, and 
only a ~ 15 % discrepancy remains between them. This clearly demonstrates that weak 
lensing tomography can break the m Vl w degeneracy very effectively. 
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Table 5. Projected lex constraints on the neutrino mass from various combinations of 
Planck and an LSST-like wide lensing survey (with one and five tomography bins) for 
two cosmological models. The eleven parameter model corresponds to that presented 
in section 4.2, while the seven parameter model assumes {N e g, a s ,wo, w a } to be fixed 
at {3,0,-1,0}. 



Model 


Cosmological probes 




11 parameters 


Planck only 


0.48 eV 


11 parameters 


Planck+Wide-1 


0.15 eV 


11 parameters 


Planck+Wide-5 


0.043 eV 


7 parameters 


Planck+Wide-1 


0.082 eV 


7 parameters 


Planck+Wide-5 


0.037 eV 



5. Conclusions 

Cosmological probes are already at an accuracy level to probe important aspects of 
particle physics such as neutrino masses and the dark energy equation of state. Current 
CMB, LSS, SNIa and BAO data can already pin the former down to < 0.5 eV [15]. 

With future CMB data and weak lensing surveys, the projected sensitivity to these 
parameters will be far better. From a combined analysis of prospective data from the 
Planck mission and a future LSST-like full-sky lensing survey, the sensitivity to m u 
is expected to reach 0.05 eV, sufficient to probe the difference between the normal and 
inverted hierarchy schemes. 

It is interesting to compare our results to those derived in [54] , which also considers 
the combination of Planck data with a future weak lensing survey. For the model closest 
to our Planck+Wide-5 eleven parameter scenario (cf. Table 5), they find a sensitivity of 
o~(J2 m v) — 0.045 eV, compared to our 0.043 eV. Systematics issues in the weak lensing 
survey have been neglected in the analysis of [54]. The same is true in our analysis in 
the sense that we have assumed the systematic effects to be controlled to the accuracy 
needed in order not to seriously degrade the parameter determination. 

A similar sensitivity to m v can also be reached by combining Planck data with 
future cluster abundance measurements [99], cr(Y^m u ) = 0.034 eV, and high-redshift 
galaxy surveys [100], a{Y,m v ) = 0.059 eV (0.5 < z < 2) and 0.043 eV (2 < z < 4). Note 
that these numbers were derived in [99, 100] for an equivalent of our seven parameter 
model. However, since both methods provide a measure of the distance-redshift relation, 
and therefore share the same principle with weak lensing tomography, it is very probable 
that these alternative probes could also offer some improvement towards the degeneracy 
between the neutrino mass and the dark energy equation of state, should the latter have 
been treated as a free parameter. This remains to be verified. 

A partial breakage of the said degeneracy is also possible by cross-correlating CMB 
and LSS measurements. The analysis of [101] finds that the combination of Planck and 
future LSS data, including cross-correlation between the two data sets, can bring the 
J2 m u sensitivity to the 0.1 eV level. This, however, is not competitive with what is 
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expected from weak lensing surveys. 

Finally, the weak lensing signal can be extracted from the CMB data itself. With 
data from the Planck satellite, the authors of [102] estimate that the neutrino mass can 
be determined to o~(^m u ) = 0.51 eV without lensing extraction, a number comparable 
to our 0.48 eV (cf. Table 5). With lensing extraction the error reduces to 0.15 eV, i.e., a 
factor of three improvement. Clearly, CMB lensing adds some of the same information 
gained in a galaxy lensing survey, except that the source redshift is locked at z ~ 1100. 
Therefore, a CMB survey with a much higher resolution is needed in order to breach 
the a(J2 m u) = 0.1 eV barrier for the neutrino mass determination. With the proposed 
Inflation Probe mission, the analysis of [102] finds that the sensitivity could reach 0.04 eV 
if lensing extraction is included, a precision similar to that attainable by large scale 
tomographic galaxy lensing surveys found in our analysis. The combination of these 
techniques could plausibly lead to even better sensitivities. 
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Appendix A. The halo model 

The halo model supposes that all matter in the universe is partitioned into distinct 
units, the halos. This assumption allows one to study the universal matter distribution 
in two steps: the distribution of matter within each halo, and the spatial distribution 
of the halos themselves. In its simplest application, the halo model proposes that the 
matter power spectrum P NL (/c, z) be composed of two distinct terms (e.g., [103]), 

P NL (k, z) = P^ halo (k, z) + P 2 - halo (A;, z), (A.l) 



where 



P 1 ~ iuio (k,z)= [ dM n(M,z)(—) KA;|M)| 2 , 



52— halo 



(k,z) = P L (k,z) 



Phalo / 

/ dM n(M, z) b{M) I — ) u(k\M) 



M N 2 



(A.2) 



, Phalo, 

describe, respectively, the correlations between two points drawn from the same halo and 
from two different halos. The former dominates on small scales, while the latter rises to 
prominence on large scales and approaches the power spectrum calculated from linear 
theory P L (k, z) as k — > 0. The average matter density phaio counts all matter clustered 
in halos. For a basic ACDM set-up, phaio is well approximated by phaio — Pm = ^mPcrit, 
where p cr i t is the present critical density. 

Three pieces of information are required to complete the model. 

(i) The mass function n(M,z), usually written as n(M,z) dM = Ph a i /M /( z/ ) du, 
specifies the number density of halos of mass M at redshift z. Here, f{v) is a 
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universal function of the peak height v = z), where <5 SC = 1.68 is the 

linear overdensity at the epoch of spherical collapse, and 

at(M,z) ee / ^P^(k,z)\W(kR)\ 2 (A.3) 

is the rms linear fluctuations filtered with a tophat function W(x) = 3/x 3 (sin:r — 
xcosx) on a scale of R = ftM/A-npm) 1 / 3 . A number of variants for f{y) exists in 
the literature. Here, we adopt the version proposed by Sheth and Tormen [104]. 

(ii) The linear bias b(v) parameterises the clustering strength of halos relative to the 
underlying dark matter, and obeys / dv f(y) b(u) = 1 by construction. 

(iii) The function u(k\M) is the Fourier transform of the matter distribution p(r\M) 
within a halo, normalised to the halo mass M. A natural choice for p(r\M) would 
be the Navarro-Frenk-White profile [105, 106]. 

To generalise the halo model to mixed dark matter cosmologies, we assume that 
only baryons and CDM cluster in halos, while the neutrinos remain an almost smooth 
component. Thus, equations (A.l) to (A.3) are evaluated for P L (k,z) = P^ +C (k,z), 
where P^ +C (k, z) is the combined baryon and CDM linear power spectrum defined in 
(3.5). Correspondingly, we take the average matter density in halos to be phaio — 
(Qb + Q c ) Pcrit- The corrected spectrum P^ c (k, z) is then added with the linear neutrino 
spectrum P^(k, z) to form the total matter power spectrum P m (k, z) as per (3.4). 
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